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ABSTRACT 

We present theoretical atmosphere, spectral, and light-curve models for extrasolar giant planets 
(EGPs) undergoing strong irradiation for which Spitzer planet/star contrast ratios or light curves have 
been published (circa June 2007). These include HD 209458b, HD 189733b, TrES-1, HD 149026b, HD 
179949b, and v And b. By comparing models with data, we find that a number of EGP atmospheres 
experience thermal inversions and have stratospheres. This is particularly true for HD 209458b, HD 
149026b, and v And b. This finding translates into qualitative changes in the planet/star contrast 
ratios at secondary eclipse and in close-in EGP orbital light curves. Moreover, the presence of atmo- 
spheric water in abundance is fully consistent with all the Spitzer data for the measured planets. For 
planets with stratospheres, water absorption features invert into emission features and mid-infrared 
fluxes can be enhanced by a factor of two. In addition, the character of near-infrared planetary spectra 
can be radically altered. We derive a correlation between the importance of such stratospheres and 
the stellar flux on the planet, suggesting that close-in EGPs bifurcate into two groups: those with 
and without stratospheres. From the finding that TrES-1 shows no signs of a stratosphere, while 
HD 209458b does, we estimate the magnitude of this stellar flux breakpoint. We find that the heat 
redistribution parameter, P„, for the family of close-in EGPs assumes values from ^0.1 to ~0.4. This 
paper provides a broad theoretical context for the future direct characterization of EGPs in tight 
orbits around their illuminating stars. 

Subject headings: stars: individual (HD 209458, HD 189733, TrES-1, HD 149026, v And, HD 179949)— 
(stars:) planetary systems — planets and satellites: general 



I 1. INTRODUCTION 

To date, more than 250 extrasolar planets have been 
discovered and more than 29 of them are transiting their 
primary star"^. One transiting planet is a "Neptune" (GJ 
436b), but the rest are giant planets with an impressively 
[ wide range of masses and radii that speak to the hetero- 
geneity of the family of close-in EGPs (Extrasolar Giant 
Planets). Table [1] lists these transiting EGPs and some 
of their relevant properties, along with many of the ref- 
erences to the observational and discovery papers from 
which these data were taken. Table [2] lists useful data for 
the corresponding primary stars, including their masses, 
luminosities, radii, and approximate ages. Both tables 
are in order of increasing planetary semi-major axis and, 
considering the pace of the field, both should be consid- 
ered provisional. Not shown are the eccentricities, which 
are generally small, but which for HAT-P-2b, GJ 436b, 
and X0-3b are -0.507, -0.14, and ~0.22, respectively. 
For these three close-in EGPs, significant tidal heating 
and, perhaps, forcing by an unseen companion are im- 
phed. 

Radial-velocity measurements for a non-transiting 
EGP provide a lower limit to its mass, but little else. 
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However, the transiting EGPs yield radii as well, and re- 
solve the sini ambiguity to reveal the planets' masses. 
These data provide physical constraints with which de- 
tailed evolutionary and structural models that incorpo- 
rate irradiation and migration can be tested (see, e.g.. 
Burrows et al. 2007a; Guillot et al. 2006). With su- 
perb photometric accuracy, the wavelength-dependence 
of the transit radii can in principle provide a measure 
of a planet's atmospheric composition (Fortney et al. 
2003). In this way, sodium has been detected in HD 
209458b (Charbonneau et al. 2002) and water has been 
identified in both HD 189733b (Tinetti et al. 2007, but 
see Ehrenreich et al. 2007) and HD 209458b (Barman 
2007). Moreover, high- precision optical photometry has 
constrained (perhaps, measured) the geometric albedo of 
the close-in EGP HD 209458b (Rowe et al. 2006, 2007). 
HD 209458b's optical albedo is very low (-3.8±4.5%), 
in keeping with the predictions of Sudarsky et al. (2000) 
when the alkali metals, and not clouds, dominate ab- 
sorption in the atmosphere and Rayleigh scattering dom- 
inates scattering. 

Nevertheless, using current technology, transit mea- 
surements have limited utility in characterizing the atmo- 
spheres and compositions of these planets. Astronomers 
require more direct detections of the planet's spectrum to 
probe its chemistry and atmospheric properties. This is 
in the tradition of remote sensing in the solar system and 
of the Earth. Until recently, it had been thought that the 
light from an extrasolar planet had to be separated from 
under the glare of its parent star using high-angular reso- 
lution, extremely high-contrast imaging. This is still the 
case in the optical for the cool "wide-separation" EGPs 
(Burrows et al. 2004; Sudarsky et al. 2005; Burrows 
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2005) and terrestrial planets, for which the planet-star 
contrast ratio is ~10^^ to ^10^^°, but such performance 
has not yet been demonstrated. 

However, for the hot close-in EGPs, the planet-star 
contrast ratios in the mid-infrared are much more fa- 
vorable (Burrows, Sudarsky, & Hubeny 2003,2004), of- 
times exceeding 10~^. This capability has led to a break- 
through in the study of exoplanets. With the infrared 
space telescope Spitzer (Werner & Fanson 1995), using 
its IRAC and MIPS cameras and the IRS spectrome- 
ter, one can now measure the summed light of the planet 
and the star in and out of secondary eclipse, and thereby, 
from the difference, determine the planet's spectrum at 
superior conjunction. Moreover, for a subset of the clos- 
est EGPs it is possible to use Spitzer to measure their 
flux variations with planetary phase between transit and 
secondary eclipse. Hence, for the close-in EGPs in the 
near- to mid-infrared, and without the need to separately 
image planet and star, the direct detection of planetary 
atmospheres via low-resolution spectroscopy and preci- 
sion infrared (IR) photometry is now a reality. 

The secondary eclipse fluxes have now been measured 
for five transiting EGPs (HD 189733b, TrES-1, HD 
209458b, HD 149026b, GJ 436b), but not yet in aU 
Spitzer bands. In addition, using the IRS spectrome- 
ter, spectra between '^7.5 ^m and ^15 /im of two tran- 
siting EGPs at secondary eclipse have been obtained 
(HD 189733b [GriUmair et al. 2007] and HD 209458b 
[Richardson et al. 2007]). Though at very low resolu- 
tion, these are the first measured spectra of any extraso- 
lar planet. Furthermore, light curves have been measured 
for three EGPs {v And b [Harrington et al. 2006: at 24 
^m], HD 179949b [Cowan et al. 2007: at 8 /xm], and HD 
189733b [Knutson et al. 2007b: at 8 /zm]). Only one 
of these (HD 189733b) is transiting and has an absolute 
calibration. For none of the latter three are there ex- 
tant light-curve measurements for more than one Spitzer 
band; for some of these EGPs only upper limits in a few 
of the other bands have been determined. Table [3] sum- 
marizes all the direct detection data for the EGP family 
obtained to date (circa June 1, 2007), along with associ- 
ated references, comments, and table notes. Clearly, in 
the next year or two we can expect a great deal more sec- 
ondary eclipse and light-curve data in the various Spitzer 
IRAC and MIPS bands. However, there has already been 
significant progress in measuring EGP atmospheres. 

This paper is a continuation of our series of interpre- 
tative studies (e.g.. Burrows et al. 2005, 2006, 2007b) of 
the direct measurements of close-in EGPs. Here, we an- 
alyze the secondary eclipse and light-curve data summa- 
rized in TableOfor the EGPs HD 189733b, HD 209456b, 
TrES-1, HD 149026b, HD 179949b, and v And b and 
make theoretical predictions in support of future Spitzer 
planet measurement campaigns. Importantly, by fitting 
the current data we extract physical information concern- 
ing the atmospheres, compositions, and thermal profiles 
of these first six objects listed in Table [H We have ex- 
plored the dependence of the spectra and light curves on 
the heat redistribution factor P„ (Burrows et al. 2006), 
on atmospheric metallicity, and on the possible presence 
of a stratospheric absorber. The recent analysis by Bur- 
rows et al. (2007b) of the IRAC data of Knutson et al. 
(2007c) indicates that HD 209458b boasts a thermal in- 
version that radically alters the Spitzer fluxes and their 



interpretation. In fact, in Burrows et al. (2007b), we 
speculate that thermal inversions and stratospheres may 
play a role in the atmospheres of many close-in EGPs 
and are a new feature in the study of transiting plan- 
ets. A similar conclusion was reached by Fortney et al. 
(2006), specifically in the context of HD 149026b. 

We find that we can fit all the secondary eclipse and 
light-curve data, except for the nightside flux of HD 
189733b and its day/night contrast. While we can fit 
its dayside secondary eclipse flux, we suspect that HD 
189733b will require a more sophisticated day/night re- 
distribution model than we now employ (32]) ■ We note 
that Knutson et al. (2007b) conclude that the dimmest 
and brightest spots on HD 189733b reside on the same 
hemisphere and that the dimmest spot is shifted from 
the anti-stellar point by as much as ~30°. Our current 
light curve models are symmetric about the peak. 

We find that the degree of longitudinal heat redistri- 
bution (P„) may vary from planet to planet, hinting at a 
variety of meteorological conditions and day/night con- 
trasts within the family of close-in EGPs. Moreover, as 
also concluded in Burrows et al. (2005), we can not ob- 
tain good fits at secondary eclipse without the presence 
of water in abundance in the atmospheres of these irra- 
diated EGPs. This is particularly true for TrES-1 and 
HD 209458b. Though we will not dwell in this paper 
on metallicity, we find that the metallicity dependence 
of the secondary-eclipse fluxes is not strong, but that 
it is in principle measurable, and that the metallicity 
dependence of the variation of the planetary flux with 
phase is only modest. Importantly, we also conclude that 
upper-atmosphere absorption in the optical by an as-yet 
unknown molecule, and the concommitent thermal in- 
versions and stratospheres, provide better fits to some of 
the data. 

In f}2l we describe our numerical techniques, the new 
heat redistribution model, and how we generate strato- 
spheres. This section is supported with Appendices iRl 
i jBl i }Cl and IJO] in which we provide details concern- 
ing the heat redistribution model and derive some ana- 
lytic formulae concerning atmospheric physics with day- 
night coupling. In particular, in fjDl we address the en- 
hancement at secondary eclipse in the integrated plan- 
etary flux at Earth over and above what would be ex- 
pected if the planet emitted isotropically. For a radar 
antenna, this would be its "gain" factor. In ij3l we dis- 
cuss the derived temperature-pressure profiles on the day 
and the night sides for all six EGPs highlighted in this 
investigation. We show that the t = 2/3 decouphng 
layers for the Spitzer IRAC and MIPS (24-^m) band 
fluxes are all above the isothermal region of an irradi- 
ated EGP's atmosphere and, hence, that Spitzer does 
not probe these deeper regions. For the same reasons, 
we find that the presence of a thermal inversion at al- 
titude and of a stratosphere in some EGP atmospheres 
can significantly alter these Spitzer fluxes and their rel- 
ative strengths. In S}3l we also provide representative 
planet spectra to demonstrate that most of the planet's 
flux emerges at shorter wavelengths than are accessible 
to Spitzer, and, hence, that Spitzer probes only a small 
tail of the emergent flux distribution. This may be of 
relevance when JWST is available to follow up on the 
Spitzer EGP data and, even earlier, as the JWST exo- 
planet campaign is being designed. 
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In ijU we present the best-fit planet-star contrast ra- 
tios at secondary eclipse for four of the transiting EGPs 
for which these have been measured (all but GJ 436b, for 
which see Deming et al. 2007 and Demory et al. 2007), as 
well as various comparison models to gauge a few param- 
eter dependences. Then, in ^we match our theoretical 
phase light curves with the three measured light curves 
and derive approximate planetary parameters. The pa- 
per is brought to a close in [JS] with a synopsis of our 
results and a general discussion of the issues raised. 

2. NUMERICAL TECHNIQUES 

Our model atmospheres are computed using the up- 
dated code CoolTLUSTY, described in Sudarsky, Bur- 
rows, & Hubeny (2003), Hubeny, Burrows, & Sudarsky 
(2003), and Burrows, Sudarsky, & Hubeny (2006), which 
is a variant of the universal spectrum/ atmosphere code 
TLUSTY (Hubeny 1988; Hubeny & Lanz 1995). The 
molecular and atomic opacities are taken from Sharp & 
Burrows (2007) and the chemical abundances, which in- 
clude condensate rainout, are derived using the thermo- 
chemical model of Burrows & Sharp (1999), updated as 
described in Sharp & Burrows (2007) and in Burrows et 
al. (2001). 

To handle convection, we use standard mixing-length 
theory, with a mixing length equal to the pressure scale 
height. The stellar irradiation boundary condition is nu- 
merically challenging, and has not been done properly by 
some workers in the past. To ensure an accurate numer- 
ical solution with a non-zero incoming specific intensity, 
we use the formalism discussed in Hubeny, Burrows, & 
Sudarsky (2003). The stellar spectral models are taken 
from Kurucz (1994) for the six stars listed in Table[2]that 
are the primaries of the EGPs upon which we focus in 
this paper. The day and night sides are approached dif- 
ferently, with the nightside, quite naturally, experiencing 
no incident flux, but receiving heat from the irradiated 
dayside using the new algorithm described in Appendix 
i Rl An important additional feature of our new heat re- 
distribution formalism is the capacity to match both the 
entropy and the gravity at the base of both the day- and 
the night-side atmospheres. Since the inner convective 
zone, which constitutes most of the planet, is isentropic, 
this is the physically correct procedure. For a dayside 
calculation, we can assume a given interior flux effec- 
tive temperature. Tint (a standard number could be 75 
K). For a given gravity and irradiation regime, this leads 
to an atmosphere solution on the dayside. This solu- 
tion incorporates an entropy in the convective zone. For 
the nightside atmosphere, we can adjust the nightside 
Tint until the entropy in its convective interior matches 
that found in the dayside convective zone. One prod- 
uct of this procedure is a connection between dayside 
and nightside Tint that has a bearing on overall planet 
cooling and shrinkage (Burrows et al. 2007a). However, 
since for a given measured planet radius, this mapping 
does not have a significant effect on the close-in planet's 
spectrum, we do this here only approximately and leave 
to a later paper a general discussion of this topic. 

The simple parametrization we use to simulate the 
effects of an extra stratospheric absorber entails plac- 
ing an absorber with constant opacity, Kq, in the opti- 
cal frequency range (z^Ojf^i) = (3 x lO^'' Hz, 7 x 10^^ 
Hz) and high up at altitude, where the pressure is be- 



low a prescribed value, generally take to be 0.03 bars. 
Hence, Hc is the most important parameter to be ad- 
justed. We could easily introduce a specified frequency 
and/or depth dependence, but this would add free pa- 
rameters which we feel are not justified at this stage. 
We have also generated models in which TiO and VO 
are allowed to assume their chemical-equilibrium upper- 
atmosphere abundances (Sharp & Burrows 2007), uncor- 
rected for the cold-trap effect (f^Hl Burrows et al. 2007b), 
and to generate a stratosphere. These TiO/VO models 
produce qualitatively the same effects as do our ad hoc 
models. However, in this paper we prefer the flexibility 
of the Ko prescription. 

In Burrows et al. (2006), once the day- and night- 
side atmospheres were calculated, we used a 2D radia- 
tive transfer code to determine the integrated emissions 
"at infinity" at a given viewing angle from the planet- 
star axis for the day- and night-side hemispheres. These 
numbers were then transformed into a light curve as a 
function of wavelength and planetary phase angle using 
the methods described in Sudarsky et al. (2000, 2005). 
However, we have found that since the detailed shape of 
the light curve connecting the day- and the night-side 
fluxes is likely to be only poorly constrained for the fore- 
seeable future and since our model is symmetric about 
secondary eclipse, it is inappropriate to invest a dispro- 
portionate amount of effort in performing expensive 2D 
transfer calculations. Rather, we invest our efforts in ob- 
taining state-of-the-art day- and the night-side fluxes at 
secondary and primary eclipse and then connect them 
with a simple, though well-motivated, curve. Therefore, 
our light curve model is: 

^ D+N D-N 

C= \ cosasmi, (1) 

2 2 ^ ^ 

where C is the planet/star flux ratio, D is the dayside 
flux (see Fig. [4]), is the nightside flux, a is the phase 
angle, and i is the inclination angle (~ 90° if in transit). 
This is the form adopted in Cowan et al. (2007). Using 
this simpler approach, one does not imply more precision 
than is warranted at this preliminary stage of inquiry. 

We have revisited the question of what type of av- 
eraging of the incoming radiation from the star over 
the surface of the planet is best suited to describe the 
planetary spectrum close to the secondary eclipse. We 
had demonstrated earlier (Sudarsky et al. 2005) that 
detailed 2D phase-dependent spectra averaged over the 
phase are equal, within a few percent, to the spectrum 
computed for a representative model atmosphere that 
is constructed assuming that the incoming flux is dis- 
tributed evenly over the surface of the dayside. In the 
usual terminology, this corresponds to the flux distribu- 
tion factor / = 1/2 (Burrows et al. 2000)^^. 

However, the spectrum of a planet observed close to 
the secondary eclipse should be biased toward a higher 
flux than that obtained using the / = 1/2 model. This 
is because the hottest part of the planet, the substellar 
point, is seen as emerging from the planet perpendicu- 
larly to the surface, and, thus, with the lowest amount 
of limb darkening. Therefore, the contribution of the 

■* We allow the planet to be irradiated by the full flux received 
from the star; it is only deeper in the atmosphere where the energy 
is transported to the nightside. See Appendix lAl 



4 



hottest part of the planet is maximized. To study this 
effect, we have computed a series of model atmospheres 
corresponding to a number of distances from the sub- 
stellar point, and have integrated the individual contri- 
butions to get the flux received by an external observer 
at a phase close to superior conjunction. It turns out 
that the flux is very close to that computed for / = 2/3, 
which is the value we subsequently use in all simulations 
presented here. There is a simple analytic argument why 
/ should be approximately equal to 2 /3 which we present 
in Appendix fjD] 

The formalism for D and N is the best we have fielded 
to date. Nevertheless, a 3D general circulation model 
(GCM) that incorporates state-of-the-art opacities, com- 
positions, and radiative transfer will be needed to prop- 
erly address day-night heat redistribution, the vortical 
and zonal mass motions, and the positions of the hot 
and cold spots. Such a model is not yet within reach, 
but there have been preliminary attempts to treat this 
physics (Showman & Guillot 2002; Cho et al. 2003; 
Menou et al. 2003; WiUiams et al. 2006; Cooper & Show- 
man 2005; Lunine & Lorenz 2002). The proper GCM 
physics remains the major uncertainty in current plane- 
tary secondary eclipse and light-curve modeling. 

3. TEMPERATURE-PRESSURE PROFILES 

We have used the techniques outlined in 321 and in Ap- 
pendix iRland the data in Tables [T] and [2] to create mod- 
els of six of the close-in EGPs in Tabic [3] for which there 
are Spitzer secondary eclipse or phase light curve data. 
These planets are HD 209458b, HD 189733b, TrES-1, 
HD 179949b, HD 149026b, and v And b. The product 
of our investigation is an extensive collection of atmo- 
sphere models, with associated spectra, for many com- 
binations of planet, P„, metallicity, and values of Pq/Pi 
We have, however, settled on presenting in this pa- 
per only the central and essential results for each planet, 
in the knowledge that the data are not yet exquisitely 
constraining. 

We focus on models with solar-metallicity (Asplund, 
Grevesse, & Sauval 2006) opacities, (Pq, -Pi) = (0.05, 0.5) 
bars, and an interior flux Tint of 75 K. These are our 
baseline model parameters. For a given measured planet 
radius, the dependence of the models on Tint is extremely 
weak. We flnd that the speciflc pressure range (Po,Pi) 
in which most of the heat carried from the dayside to the 
nightside is conveyed plays a role in the planet-star flux 
ratios, but a subtle one. Therefore, in lieu of a compre- 
hensive, and credible, 3D climate model, we prefer not 
to claim too much concerning the details of atmospheric 
circulation and heat redistribution. We also explore the 
effects of a stratospheric absorber with an optical opacity 
of Ke (see also Burrows et al. 2007b). We find that such 
models will be most important for close-in EGPs with 
the greatest stellar insolation and guided by this princi- 
ple, particularly relevant for HD 209458b, HD 149026b, 
and V And b, we explore the consequences. 

Figure [1] portrays in six panels the temperature- 
pressure (T/P) profiles of a representative collection of 
dayside and nightside models of the six close-in EGPs 
of this study. For the dayside, the different curves cor- 
respond to different values of P„ from 0.0 (no redistri- 
bution) to 0.5 (full redistribution) and to models with 
and without stratospheric absorbers. For the nightside. 



P„ ranges from 0.1 to 0.5. For all models, the radiative- 
convective boundaries are identified and are quite deep 
(on the far right of each panel). When Kc ^ 0, the T/P 
profiles show distinct thermal inversions. 

There are quite a few generic features in evidence on 
these panels. The first is that the atmospheres are never 
isothermal. Since the opacities in the optical, where 
most of the stellar irradiation occurs, and in the infrared, 
where most of the reradiation occurs, are very different, 
a quasi- isothermal inner region interior to '^l bar is al- 
ways bounded by an outer region in which the tempera- 
ture decreases (Hubeny, Burrows, & Sudarsky 2003). As 
Fig. [1] indicates, the magnitude of the temperature de- 
crease from the plateau to the ~10~^ bar level for dayside 
models without stratospheres is ■^lOOO K. With a strato- 
sphere, the outward increase from a pressure of ~0.1 bars 
can be correspondingly large. For our nightside models, 
the monotonic decrease is ~500— 1000 K. Models with 
temperature inversions due to a strong absorber at alti- 
tude clearly stand out in the panels of Fig. [T] and may 
result from the presence of a trace species, TiO/VO, or 
a non-equilibrium species (Burrows et al. 2007b). The 
possible effect of such upper-atmosphere absorbers on 
the T/P profiles and the resultant dayside spectra are 
exciting new features of the emerging theory of irradi- 
ated EGPs. 

The discussion above is made more germane when we 
note that the decoupling surfaces for the IRAC and MIPS 
(24-/xm) channels, the effective photospheres where t\ ^ 
2/3, are all in the outer zone. Figure [1] indicates their po- 
sitions for the dayside P„ — 0.3 model of TrES-1. They 
are at similar pressures for all other models. The photo- 
spheres for shorter wavelengths not accessible to Spitzer 
are deeper in. Figure [2] portrays these "formation," 
"brightness," or photospheric temperatures as a function 
of wavelength for three models of TrES-1 with different 
values of P„, and for both the dayside and nightside, and 
illustrates this fact. The approximate wavelength inter- 
vals of the Spitzer bands are superposed. Though the 
IRAC 1 flux can decouple at interesting depths, the pho- 
tospheres for the Y, J, H, and K bands are generally 
deeper. The photospheres in the far-IR beyond ~10 /im 
are high up at altitude and we repeat that Spitzer pho- 
tometry does not probe the isothermal region so charac- 
teristic of theoretical close-in EGP atmospheres. More- 
over, since the Spitzer observations are probing the outer 
regions of the atmosphere most affected by stellar irradia- 
tion and which can have inversions, the treatment of the 
outer boundary condition due to incoming stellar flux 
must be accurate. Slight errors or uncertainties in the 
outer boundary condition of the transfer solution, or in 
the upper-atmosphere opacities, can translate into sig- 
nificant errors in the predicted Spitzer fiuxes. This is 
particularly true longward of ~10 /im. As a result, mea- 
sured fluxes in both the near-IR and mid-IR are useful 
diagnostics of upper-atmosphere absorbers and thermal 
inversions (Hubeny, Burrows, & Sudarsky 2003; Burrows 
et al. 2006; Burrows et al. 2007b; Knutson et al. 2007c). 
All these caveats and points must be borne in mind when 
interpreting the Spitzer EGP data. 

As a prelude to our discussions in S|4] and fj5] of the 
Spitzer planet-star flux ratios at secondary eclipse and 
during an orbital traverse, and to emphasize the fact that 
Spitzer does not comprehensively probe the irradiated 
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planet's atmosphere, we plot in Fig. [3]theoretical dayside 
spectra (AFa versus logio(A)) for three models of TrES-1 
at zero phase angle (superior conjunction). Superposed 
on the plot are the positions of the near-IR, IRAC, and 
MIPS bands. Such a figure allows one to determine at 
a glance the wavelengths at which most of the flux is 
radiated (at least, theoretically). As Fig. [3] suggests, 
most of the planet's flux emerges in the near-IR, not in 
the IRAC or MIPS channels. In fact, depending on the 
planet, no more than one fifth to one third of the planet's 
flux comes out longward of ^3.6 fim (IRAC 1), while no 
more than one twentieth to one tenth emerges longward 
of ~6.5 /im, the "left" edge of the IRAC 4 channel. Since 
much of the best EGP data have been derived in IRAC 
channel 4, one must acknowledge that they may represent 
very little of the total planetary emissions. 

Finally, we call the reader's attention to the slight 
bumps (on the nightside) and depressions (on the day- 
side) between 0.05 and 1.0 bars in the T/P profiles de- 
picted in Fig. [TJ This region is near where we imposed 
heat redistribution using the formalism described in iRl 
The actual shapes of these profiles are determined by this 
mathematical procedure and other algorithms will pro- 
duce different local thermal profiles. Note that with this 
formalism it is possible at the higher P„s (> 0.35) for the 
nightside to be hotter than the dayside at the same pres- 
sure levels in the redistribution region. While this may 
seem at odds with thermodynamics, what is essential is 
that energy is conserved and is redistributed at optical 
depths that are not either too low or too high. If the 
former, the absorbed stellar heat would be reradiated 
before it can be carried to the nightside. If the latter, 
then the stellar radiation can not penetrate to the con- 
veyor belt. For our default choice of Pq and Pi, TRossiand 
is generally between ~0.3 and ~6. These depths are not 
unreasonable, but our redistribution algorithm is clearly 
only a stopgap until a better GCM can be developed and 
justified. 

4. PLANET-STAR FLUX RATIOS-COMPARISON WITH 

DATA 

We discuss below and in turn model fits for each 
transiting EGP at secondary eclipse. However, first we 
present our results collectively and in summary fashion. 
Figure [5] in four panels portrays for the four transiting 
EGPs the correspondence between the secondary eclipse 
data and representative models of the planet-star flux 
ratio as a function of wavelength from 1.5 /im to 30 fj,m. 
This figure summarizes our major results. The models 
are for values of P„ of 0.1, 0.3, and 0.5 and various values 
of He- The data include l-cr error bars and can be found 
in Table [31 As Fig. [3] indicates, we can fit all the pub- 
lished data. The P„ dependence for both stratospheric 
models and models without inversions is strongest in the 
K band and in IRAC 1. In fact, in the near-IR, models 
with inversions depend very strongly on P„. Fig. [5] for 
HD 209458b in the near-IR indicates this most clearly. 
This finding implies that measurements at these shorter 
IR wavelengths are good diagnostics of P„, particularly 
if inversions are present. 

Importantly, including a non-zero Hc and generating a 
stratosphere results in a pronounced enhancement long- 
ward of IRAC 1, particularly in IRAC 2 and 3, but 
also at MIPS/24 /im and at the 16-/im peak- up point 



of Spitzer/IRS. Hence, fluxes at the longer IR wave- 
lengths might be good diagnostics of thermal inversions. 
The models in Fig. H for HD 209458b and HD 149026b 
demonstrate this feature best. 

No attempt has been made to achieve refined fits, but 
the correspondence between theory and measurement, 
while not perfect, is rather good for all the planets. More- 
over, different EGPs seem to call for different values of 
P„ and Ke, and, hence, perhaps, different climates, de- 
grees of heat redistribution, compositions, and upper- 
atmospheric physics. The light-curve analyses in [J5] also 
suggest this. A goal is to relate these measured differ- 
ences with the physical properties of the star and planet 
and these infrared secondary eclipse data allow us to be- 
gin this program in earnest. 

We note that comparisons between model and data 
must actually be made after the band-averaged flux- 
density ratios of the detected electrons are calculated. 
Performing this calculation slightly mutes the predicted 
variation from channel to channel in the IRAC regime. 
This is particularly true when comparing IRAC 1 and 
IRAC 2, even if a pronounced spectral bump at and near 
~3.6 /im obtains, as it does for models with modest or 
no thermal inversion. However, to avoid the resultant 
clutter and confusion, we do not plot these bandpass- 
averaged predictions on Figs. [J] and [5l We now turn 
to case-by-case discussions of the secondary eclipse mea- 
surements and models. 

4.1. HD 209458b 

The first transiting EGP discovered was HD 209458b 
(Henry et al. 2000; Charbonneau et al. 2000) and it has 
since been intensively studied. The direct-detection data 
of relevance to this paper are summarized in Table[3] The 
most relevant data are the geometric albedo constraints 
in the optical from MOST (Rowe et al. 2006,2007), a 
i^-band upper limit using IRTF/SpeX from Richard- 
son, Deming, & Seager (2003), a MIPS/24-/^m photo- 
metric point from Deming et al. (2005) (and its possi- 
ble update) , a low-resolution Spitzer/IRS spectrum from 
Richardson et al. (2007), and, importantly, photomet- 
ric points in IRAC channels 1 through 4 from Knutson 
et al. (2007c). These data collectively provide useful 
information on the atmosphere of HD 209458b. 

Motivated by the recent data of Knutson et al. (2007c), 
Burrows et al. (2007b) provide partial theoretical ex- 
planations for HD 209458b's atmosphere. Much of the 
discussion in the current paper concerning HD 209458b 
is taken from Burrows et al. (2007b), so we refer the 
reader to both the Burrows et al. (2007b) and Knutson et 
al. (2007c) papers for details. However, here we expand 
upon the discussion in those works where it is necessary 
to put the HD 209458b findings in the broader context 
of the EGPs listed in Table [3l The major conclusion of 
Burrows et al. (2007b) is that the atmosphere of HD 
209458b has a thermal inversion and a stratosphere, cre- 
ated by the absorption of optical stellar flux by a strong 
absorber at altitude, whose origin is currently unknown. 
This converts absorption features into emission features, 
while still being consistent with the presence of water in 
abundance. 

All relevant data, save the albedo constraint in the 
optical, are displayed in the upper-left panel of Fig. 21 
Figure [5] includes the Knutson et al. (2007c) IRAC 1 
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point and the Richardson, Deming, & Seager (2003) up- 
per limit in K and focuses on the near-IR. Also provided 
on both figures are models for P„ = 0.1, 0.3, and 0.5, 
without and with an extra stratospheric absorber. The 
latter is implemented using the formalism in outlined in 
fj2]and a Ko of 0.1 cm^/g. 

Figure S] shows that the low upper limit of Richard- 
son, Deming, & Seager (2003) in the K band that 
was problematic in the old default theory (Burrows, 
Hubeny, & Sudarsky 2005; Fortney et al. 2005; Barman, 
Hauschildt, & AUard 2005; Seager et al. 2005; Burrows 
et al. 2006) is consistent with the models with an extra 
upper-atmosphere absorber in the optical, particularly 
for higher values of P„. This is more clearly seen in Fig. 
[5l Moreover, the theoretical peak near the IRAC 1 chan- 
nel (~3.6 /im) in the old model without an inversion is 
reversed with the extra absorber into a deficit that fits 
the Knutson et al. (2007c) point. The theory without 
an extra absorber at altitude predicts that the planet- 
star flux ratio in the IRAC 2 channel should be lower 
than the corresponding ratio at IRAC 1. However, with 
the extra absorber the relative strengths in these bands 
are reversed, as are the Knutson et al. (2007c) points. 
This reversal is a direct signature of a thermal inversion 
in the low-pressure regions of the atmosphere, and an 
indirect signature of the placement of the heat redistri- 
bution band (see Appendix fjA|. The top- left panel of 
Fig. [1] depicts the corresponding temperature-pressure 
profiles and the thermal inversion at low pressures intro- 
duced by the presence of an extra absorber in the optical 
that is indicated by the data. 

As Fig. [4] also demonstrates, there is a significant dif- 
ference in the IRAC planet-star flux ratios between the 
old default model without an inversion and the new mod- 
els with a stratospheric absorber, and that the models 
with a stratosphere fit the IRAC channel 1, 2 and 4 flux 
points much better. However, the height of the IRAC 3 
point near 5.8 fiui is not easily fit, while maintaining the 
good fits at the other IRAC wavelengths and consistency 
with the if-band limit. Theoretically, the positions of the 
IRAC 3 and IRAC 4 photospheres should be close to one 
another, so this discrepancy is surprising. Nevertheless, 
the IRAC 2, 3, and 4 data together constitute a peak, 
whereas in the default theory an absorption trough was 
expected. 

The 24-/xm MIPS point obtained by Deming et al. 
(2005) is lower than the prediction of our best-fit model. 
However, the fiux at this point is being reevaluated and 
may be closer to ~0.0033±0.0003 (D. Deming, private 
communication). If the new number supercedes the old 
published value, then our best-fit model(s) with inver- 
sions fit at this mid-IR point as well (Fig. 0]). Higher 
planetary fluxes longward of ^10 /im are generic features 
of stratospheric inversions. 

The 1-cr optical albedo limit from Rowe et al. (2007) 
is a very low 8.0%. For comparison, the geometric albe- 
dos of Jupiter and Saturn are ~40%. However, such 
a low number was predicted due to the prominence in 
the optical of broadband absorptions by the alkali metals 
sodium and potassium in the hot atmospheres of irradi- 
ated EGPs (Sudarsky et al. 2000; their "Class IV"). The 
associated planet-star flux ratios are ~10~^ — 10"^. This 
low albedo is consistent with the identification of sodium 
in the atmosphere of HD 209458b using HST/STIS tran- 



sit spectroscopy (Charbonneau et al. 2002). Both these 
datasets suggest that any clouds that might reside in the 
atmosphere of HD 209458b are thin. A thick cloud layer 
would reflect light efficiently, leading to a high albedo. If 
the extra stratospheric absorber is in the gas phase, and 
there is no cloud, then our new thermal inversion mod- 
els are easily consistent with the low albedo derived by 
Rowe et al. (2006,2007). If the extra absorber is a cloud, 
the cloud particles must have a low scattering albedo and 
can not be very reflecting. This rules out pure forsterite, 
enstatite, and iron clouds. 

The IRS data are noisy, but their flattish shape is con- 
sistent with our model(s) with thermal inversions and a 
stratosphere. Richardson et al. (2007) suggest that there 
is evidence in the IRS data for two spectral features: one 
near 7.78 /im and one near 9.67 /im. However, we think 
the data are too noisy to draw this conclusion. Richard- 
son et al. (2007) also suggest that the flatness and ex- 
tension of their data to shorter wavelengths implies the 
near absence of water, since previous theoretical mod- 
els predicted a spectral trough between ~4 /^m and ^8 
/im. However, if there is an outer thermal inversion, as 
we here and in Burrows et al (2007b) argue is the case 
for HD 209458b, a trough is flipped into a peak for the 
same water abundance. This renders moot the use of the 
spectral slope at the edge of the IRS spectrum to deter- 
mine the presence or absence of water. One of our major 
conclusions, implicit in Figs. [3] and [5l is that water is 
not depleted at all in the atmosphere of HD 209458b. 

The recent controversies surrounding such an interpre- 
tation occasion the following remarks. One thing to bear 
in mind concerning the use of these IRS spectra to infer 
compositions is that they are very low resolution. The 
use of classical astronomical spectroscopy to identify con- 
stituents stems from the ability at much higher resolution 
to see characteristic features at precise wavelengths and 
patterns of absorption or emission lines to high accuracy. 
This allows one to make element and molecule identifi- 
cations in a narrow wavelength range without a global 
view across the whole spectrum. However, at the low 
resolution of the IRS, no individual water features are 
accessible. There are water features near '^10 //m, but 
they would require a A/AA is excess of ~2000 to identify. 
Otherwise, all one sees is the collective effect of millions 
of lines and the resulting pseudo-continuum (the band 
structure). Clearly, when the data are low- resolution, 
a global photometric and spectral fit is necessary to ad- 
dress the issue of composition. The signature of water's 
presence comes from the goodness of the global fit across 
the entire spectrum from the optical to the mid-IR. The 
good fit we obtain in Fig. [H lead us to conclude that 
the IRS, IRAC 4, and MIPS data for HD 209458b are 
consistent with the presence of water in abundance. 

Note that if the T/P profile were entirely flat (but 
see Fig. [ij, whatever the opacity and molecular abun- 
dances the emergent spectrum would be a perfect black 
body and would give no hint concerning composition. 
Transit spectra would then be our only reliable means 
of determining atmospheric composition (Fortney et al. 
2003; Barman 2007; Tinetti et al. 2007; Ehrenreich et 
al. 2007). However, as we have argued, there is every 
indication that the T/P profiles of strongly irradiated 
EGPs are not flat (Fig. [1]). As a result, spectral mea- 
surements of irradiated EGPs can be usefully diagnostic 
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of both composition and non-trivial thermal profiles. 

4.2. HD 189733b 

Models for HD 189733b at secondary eclipse, with and 
without an extra upper-atmosphere absorber, are por- 
trayed in the upper right-hand panel of Fig. [31 They 
include Ko — 0.0 cm^/g models with P„ = 0.1, 0.3, and 
0.5, and one Kg = 0.04 cm^/g model with P„ = 0.3. 
The IRAC 4 data at 8 fim from Knutson et al. (2007b) 
[brown], the IRS peak-up point at 16 /im obtained by 
Deming et al. (2006) [gray], and the IRS spectrum 
between ~7.5 fj,m and ~13.5 /Ltm from Grillmair et al. 
(2007) [gold] are superposed on the figure. Though data 
in the other IRAC channels and at 24 /im have been 
taken and reduced, they have yet to be published. 

As this panel indicates, the IRAC 4 point can be fit 
by models that include most values of P„, with a very 
slight preference for lower values from 0.1 to 0.3. The 
IRS data are not well-calibrated, but evince the slight 
turndown at the shorter wavelengths characteristic of at- 
mospheres with weak or no stratospheric absorber. This 
turndown is in contrast with the behavior of the Richard- 
son et al. (2007) IRS data for HD 209458b, and reinforces 
the conclusion that a thermal inversion, if present in HD 
189733b, is very slight (see the upper-right panel of Fig. 
[T]). However, the 16-/im point of Deming et al. (2006) 
is a bit higher than models with = 0, whatever the 
value of P„. This suggests that there may be some ex- 
tra heating in the upper atmosphere of HD 189733b, but 
that it is weaker than in the atmosphere of HD 209458b. 
The Kc = 0.04 cmVg model shown in the HD 189733b 
panels of Figs. [1] and H] indicates the possible magnitude 
of such stratospheric heating, if present. Note that to fit 
the 16-/im point we require a smaller value of Kc than 
used to fit the IRAC channel data for HD 209458b (kc = 
0.1 cm^/g). This may not be surprising, since, as Table 
[T] indicates, the stellar flux at the substellar point of HD 
189733b is lower by more than a factor of two than the 
corresponding number for HD209458b. Perhaps, this in- 
dicates a systematic trend for the family of strongly irra- 
diated transiting EGPs (see Tables [1] and [3|), with planets 
with the higher values of Fp possessing stratospheres and 
atmospheres with pronounced inversions. 

Be that as it may, we can predict, using the logic em- 
ployed in Sgm that if the IRAC 1 to IRAC 2 ratio turns 
out to be greater than or close to one, any thermal in- 
version in the atmosphere of HD 189733b is either not 
pronounced, or is absent. Under these circumstances, we 
certainly would then expect the "brightness" tempera- 
ture at IRAC 1 to be demonstrably higher than that at 
IRAC 2 (see, e.g., Fig. ^. Conversely, if the IRAC 
1 to IRAC 2 ratio is much less than one (as for HD 
209458b), then a thermal inversion in the atmosphere of 
HD 189733b would be strongly suggested. The same can 
be said of the IRAC 4 to IRAC 3 ratio. If the IRAC chan- 
nel 3 planet-to-star flux ratio is higher than the Knut- 
son et al. (2007b) point at 8 /im, then a stratosphere 
would be indicated for HD 189733b. Since the irradi- 
ation regime of HD 189733b is a bit more benign than 
that of HD 209458b, and given the contrast in the short- 
wavelength behavior of the IRS data for each EGP, we 
hypothesize that HD 189733b does not boast much of a 
stratosphere. Note that our HD 189733b models all have 
water in abundance and that if the unpublished IRAC 1 , 



2, and 3 data for HD 189733b prove to decrease monoton- 
ically with wavelength shortward of IRAC 4, this would 
be fully consistent with the presence of the water band 
between ~4 /xm and ~8 fxm in absorption. 

4.3. TrES-1 

Charbonneau et al. (2005) obtained IRAC 2 (~4.5 
^m) and IRAC 4 (~8.0 fj.m) data for TrES-1 and these 
data were analyzed by Burrows et al. (2005). A major 
conclusion of that paper was that water is indeed seen in 
absorption. Our new models, depicted in the bottom-left 
panel of Fig. d] with the two IRAC data points super- 
posed, reinforce this finding. As can be seen in the fig- 
ure, the models with different values of P„ (here, all with 
Kc — 0) can not easily be distinguished using these two 
IRAC points. This fact emphasizes the need to obtain 
more Spitzer photometric data to help better constrain 
the properties of the atmosphere of TrES-1. However, it 
is clear from the significant drop in planet /star fiux ratio 
from IRAC 4 to IRAC 2 that the atmosphere of TrES-1 
is qualitatively different from that of HD 209458b. In 
particular, this behavior is a signature of a strong water 
absorption trough. There are no signatures of a thermal 
inversion in the atmosphere of TrES-1, or of water in 
emission, and the old, default models with a monotonic 
temperature profile (see left-middle panel of Fig. [T|) are 
perfectly suitable. Given this, we predict that the IRAC 
1 point (when obtained) will be slightly higher than the 
IRAC 2 point. 

4.4. HD 149026b 

As is suggested by the relative values of Fp found in Ta- 
ble [T] for HD 209458b, HD 189733b, and TrES-1 and the 
different thermal profiles inferred for their atmospheres, 
there appears to be a correlation between the charac- 
ter of a planet's atmosphere and its value of Fp, or a 
related quantity (UV insolation?). This possibility is in- 
triguing, but not yet explained. In particular, we have 
yet to identify the extra stratospheric absorber in that 
subset of close-in EGPs "clearly" manifesting thermal in- 
versions. However, 1) HD 149026b's Fp is almost twice 
that of HD 209458b, 2) it has one of the hottest atmo- 
spheres among those listed in Table [T] (see middle-left 
panel of Fig. [1]), and 3) we can not fit the IRAC 4 
data point obtained by Harrington ct al. (2007) without 
a strong temperature inversion. The latter conclusion 
agrees with that of Fortney et al. (2006), who predicted 
a mid-infrared fiux for HD 149026b near the value actu- 
ally measured by Harrington et al. (2007) by allowing 
TiO/VO to reside at low pressures for the hottest at- 
mospheres (Hubeny et al. 2003). The lower- right panel 
of Fig. m depicts three models for HD 149026b, one of 
which has Ke = 0.64 cm^/g. This is much larger than 
the Ko employed to fit HD 209458b. The stratospheric 
model shown has P„ = 0.0, which minimizes the value of 
Kc necessary to fit the lone Harrington et al. (2007) data 
point, and it is the only model among the three depicted 
in Fig. |4] that does fit. Therefore, the trend in "inver- 
sion" strength" with Fp, seen in the sequence TrES-1, HD 
189733b, and HD 209458b, continues with HD 149026b. 
Not only do the EGP atmospheres grow hotter with Fp 
(a not unexpected result), but the importance of a ther- 
mal inversion and a stratosphere in explaining the extant 
data increases with it as well. Given the current paucity 
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of data for HD 149026b, we urge that HD 149026b be 
a priority target so as to help discriminate the various 
models only partially represented on the HD 149026b 
panel of Fig. |4l We predict that the pattern of the four 
IRAC flux ratios for HD 149026b will mimic that found 
for HD 209458b, and that its flux ratios from ~10 fim to 
~30 nm will comfortably exceed those of models without 
obvious thermal inversions, perhaps by large margins. 

5. LIGHT CURVES-COMPARISON WITH DATA 

Measuring the infrared planet-star contrast ratio as a 
function of orbital phase, i.e., the planet's light curve, 
provides the best constraints on the longitudinal dis- 
tribution of planetary emissions. In principle, phase- 
dependent light curves at different wavelengths can be 
inverted to determine the "brightness" temperature and 
composition distributions over the surface of the planet, 
including its night side. Contrast ratios obtained not just 
at secondary eclipse (a = 0°), but also at other phase 
angles, help reveal and quantify zonal winds and estab- 
lish their role in redistributing stellar energy (i.e., P„) 
and matter around the planet. They can help identify 
transitions at the terminator (Guillot & Showman 2002; 
Showman & Guillot 2002), shifts in the substellar hot 
spot (Cooper & Showman 2005; WiUiams et al. 2006), 
asymmetries in the thermal distributions (Knutson et al. 
2007b), and persistent atmospheric structures. In sum, 
light curve measurements probe both atmospheric dy- 
namics and the planet's climate and are the key to the 
bona fide remote sensing of exoplanets. 

Having said this, since full light curves require many 
more pointings and much more telescope time to obtain, 
and mostly address the dimmer phases of a planet's or- 
bital traverse, obtaining them is much more difhcult than 
measuring the contrast ratio at secondary eclipse. As a 
result, to date there are only three published light curves 
for irradiated EGPs (for v And b, HD 179949b, and HD 
189733b), despite numerous observational forays. All of 
these are for only one Spitzer waveband each and none 
covers a complete orbit. There do exist recent upper lim- 
its (e.g.. Cowan et al. 2007), but these are not usefully 
constraining and we will not address them here^. 

Below, we discuss the three systems for which light 
curves, however sparse, have been obtained and try to ex- 
tract physical information by comparison with our light- 
curve models (eq. [T]). Before we do so, we note the 
following. Classic light curve studies are in the opti- 
cal and measure geometric albedos, phase functions (Su- 
darsky et al. 2005), and polarizations, i.e. they mea- 
sure reflected stellar light. Albedos and polarizations are 
significantly affected by the presence of clouds, and so 
these traditional optical campaigns focus on reflection 
by condensates or surfaces. The planet-star flux ratios 
in the optical range from ^^10^^^ for EGPs at AU dis- 
tances to -10-^-10"^ for the close-in EGPs near ~0.05 
AU (Table [T]). However, in the near- and mid-IR, the 
planet-star contrasts are around ^lO^'^ (see Figs. |4]and 
[5]). These larger numbers are why Spitzer IR measure- 
ments, rather than optical measurements, have assumed 
center stage in the direct study of EGPs. The light seen 
is not reflected stellar light, but reprocessed stellar flux, 
emitted predominantly in the near- and mid-lR (Fig. [3]) 

^ However, a few of these limits are listed in Table |3] 



at the lower temperatures (^^1000— 2000 K) of the re- 
sulting planetary atmospheres. The ratio between the 
optical and IR components, and thus the relative advan- 
tage of IR measurements, is roughly the square of the 
ratio between the orbital distance and the stellar radius, 
a number near ~10^ for most of the planets listed in 
Tabled 

5.1. V And b 

Harrington et al. (2006) have measured the phase vari- 
ation at 24 fim of the planet-star contrast for the close-in 
EGP V And b (Butler et al. 1997). Since this planet is 
not transiting, we know neither Mp nor sin(i), but only 
the combination MpSin (z) (= 0.69 Mj ). Moreover, with- 
out a transit we don't have a measurement of Rp . In fact, 
the models used to fit the five (!) data points obtained by 
Harrington et al. (2006), which are not anchored by ab- 
solute calibration, depend on P„, Kc, Rp , and sin (i) (see 
eq. [1]). In addition, all interpretations hinge upon only 
the two extreme points in the Harrington et al. (2006) 
dataset. Therefore, we have too many degrees of free- 
dom to allow us to draw strong conclusions concerning 
planetary and atmospheric parameters and must make 
do with limits and general correlations. 

Figure [5] portrays in eight panels comparisons of theo- 
retical 24-/xm phase curves with the v And b data. The 
left panels contain models with — 0, and the right 
panels contain models with Ko = 0.2 cm^/g. The models 
in the top four panels have P„ = 0.0, and those in the 
bottom four panels have P„ = 0.3. Inclinations of both 
45° and 80° (near eclipse) are employed. On each panel, 
we provide models with a wide range of planetary radii. 
Since the data have no absolute calibration, we are free 
to move the data points up and down, as long as their 
relative values are maintained, and we have done so in 
an attempt to provide on each panel the best fit to the 
overall shape and the day/night difference. The corre- 
sponding T/P profiles at a = 0° (day) and a = 180° 
(night) are displayed in the lower-left panel of Fig. [TJ 

From figures like Fig. [6l we can extract general trends 
and limits. In their discovery paper, the authors noted 
that the shift of the hot spot away from the substel- 
lar point was small. They also remarked on the large 
difference from peak to trough (~0.002). Both observa- 
tions suggested that there is not much heat redistribution 
from the dayside to the nightside and that P„ is small, 
perhaps near zero. While at this stage this conclusion 
can not be refuted, our models suggest that there is a 
broader range of possible interpretations. Importantly, 
as we have noted in 21 the presence of a stratosphere 
can enhance dayside planetary fluxes in the mid-IR, and 
certainly near 24 /im, even for modest values of P„ which 
would otherwise decrease D/N (eq. [1]). v And b expe- 
riences a substellar flux, Fp, of ^-^1.3x10^ erg cm~^ s~^, 
and this is larger than that impinging upon HD 209458b. 
Therefore, we expect that the atmosphere of v And b 
will have a thermal inversion as well, thereby enhanc- 
ing the mid-IR day/night contrasts even for values of P„ 
near 0.3®. What is more, large values of Rp are becoming 

® It is possible that sin (i) is small, and, therefore, that the grav- 
ity and Mp are large. It is also possible that a large gravity can 
shift the breakpoint between EGP atmospheres with and without 
inversions. However, we suspect that Fp, more than gravity, is the 
crucial parameter in determining this bifurcation. 
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commonplace, and we can not eliminate this possibility 
for V And b. A large value of Rp increases the amplitude 
swing from day to night. Therefore, many parameter 
combinations can fit these data. 

We summarize the lessons of Fig. [6] as follows. All 
else remaining the same, a change of P„ from 0.0 to 0.3 
results in an increase in the Rp required of ~0.3— 0.5 Rj . 
Substituting = 0.2 cm^/g for = 0.0 decreases the 
Rp necessary by ^0.4—0.5 i?j . Replacing models with i 
= 45° by those with i — 80° decreases the Rp required 
by ~0.2-0.3 i?j . Specifically, a (P„, Ke, i) = (0.0, 0, 
80°) model has a required radius of ~2.0 i?j (top-left 
panel of Fig. El), while a (P„, Ke, i) = (0.0, 0.2, 80°) 
model has a required radius of ~1.5 Rj (top-right panel 
of Fig. El). And while a (P„, Kc, i) = (0.3, 0, 80°) model 
requires an RpOf ~2.5 Rj , one with Kc — 0.2 cm^/g 
requires a radius of "only" ~1.8 i?j . Clearly, introducing 
stratospheres into the mix allows P„ to assume a range 
of non-zero values for which heat redistribution would 
not be considered small. However, the implied planetary 
radius would not be small either, though it would still be 
within the currently measured range. 

Some of this degeneracy might be broken with light- 
curve measurements at many wavelengths and with a 
more rapid cadence. Furthermore, astrometric mea- 
surements of the stellar wobble can provide sin (i) and 
the planet's mass, eliminating one important ambiguity. 
Spitzer can still be used to provide the former, while the 
latter is well within reach of ground-based astronomy. 
Finally, we would be remiss if we did not mention that 
JWST will inaugurate an era of stunning photometric 
improvement (by a factor of at least ~10^) over current 
IR platforms for the study of the light curves of both 
transiting and non-transiting EGPs. 

5.2. HD 17994% 

Cowan et al. (2007) obtained a light curve in IRAC 
channel 4 (~8 /im) of another non-transiting giant 
planet, HD 179949b (Santos, Israehan & Mayor 2004; 
Wittenmyer, Endl & Cochran 2007). It has an MpSin (i) 
of -0.95 Mj and a value of Fp of -1.32x10^ erg cm'^ 
s^^). This makes it similar to the v And system, both 
in its general properties and in the limitations on what 
can be uncovered. 

The lower right-hand panel of Fig. [1] displays the day- 
side and nightside T/ P profiles of the light-curve models 
used to model these data. Figure [7] portrays the cor- 
responding eight-panel figure comparing our theoretical 
models for various combinations of P„, Ko, sin(i), and 
Rpwith the eight data points of Cowan et al. (2007). 
Figure [7] is similar to Fig. El but the thermal inversion 
models are for — 0.08 cm^/g and the range of model 
radii are different. The data for HD 179949b are impor- 
tant, but no less ratty than those for v And b. Never- 
theless, the values of Rp required to fit HD 179949b are 
systematically lower. 

We summarize our conclusions from Fig. [71 as follows. 
A change of P„ from 0.0 to 0.3 requires Rp to increase by 
~0. 1-0.2 R] . Replacing = 0.0 cmVg by = 0.08 
decreases Rp by ~0.2— 0.4 i?j . Substituting models with 
i ^ 45° for those with i = 80° decreases Rpby -0.1-0.2 
i?j . A (P„, Ke, i) = (0.0, 0, 80°) model that fits the 
data requires a radius of —1.2 i?j (top-left panel of Fig. 



O, while a (P„, k^, i) = (0.0, 0.08, 80°) model requires a 
radius of —1.0 i?j (top-right panel of Fig. |71). Whereas a 
(P„, Ke, i) = (0.3, 0, 80°) model requires an Rp of —1.6 
i?j , one with Kc — 0.08 cm^/g requires a much smaller 
radius, —1.1 i?j . These radii are unexceptional, and 
in general the presence of a stratosphere substantially 
decreases the values required to fit these data. Therefore, 
for HD 179949b we can fit the light-curve data with a 
quite reasonable combination of parameters, though the 
various degeneracies still need to be broken. 

Finally, Fig. [Tlsuggests that the shift between the tran- 
sit ephemeris and the light curve phases is not large, so 
we don't see any obvious advection downstream of the 
hot spot. This is consistent with the interpretation by 
Harrington et al. (2007) of the v And b phase curve 
(Fig. El) and may be a feature of EGPs with strato- 
spheres and/or hot upper atmospheres. While very ten- 
tative, this suggestion is reinforced by the observation 
that there is a definite displacement of the hot spot of 
HD 189733b, which seems to have a cooler upper atmo- 
sphere (Fig. [1]). However, the perception of meaningful 
differences in the displacements of hot spots could just 
as easily be false and be a consequence of having better 
data for HD 189733b. Moreover, we don't yet have a 
good model for the origin of such differences and possi- 
ble correlations with P„. Clearly, better sampled phase 
curve data would be very useful. 

5.3. HD 189733b 

Currently, the only light curve we have for a transit- 
ing EGP was obtained in IRAC 4 at —8 /im by Knutson 
et al. (2007b)'^. Not only do we have HD 189733b's 
radius (Table (H), but these light-curve data have abso- 
lute calibrations. In addition, there is dense coverage 
over a bit more than half the orbit, from just before sec- 
ondary eclipse to just after primary transit. Knutson et 
al. (2007b) derive the longitudinal dependence of the 
surface brightness and find a hot spot shifted by 16 ± 6° 
east of the substellar point, while the coolest region is 
shifted about 30° west of the anti-stellar point. Curi- 
ously, both the hot spot and the coolest spot are in the 
same hemisphere. Nevertheless, this is the first "map" 
of the surface of an exoplanet (Burrows 2007). The au- 
thors also found an indication of a nonzero, but small, 
eccentricity with ecosw = 0.0010 ± 0.0002, where cu is 
the longitude of periastron, a transit radius at 8 /im of 
1.137±0.006 Rj (slightly smaller than the optical radius), 
a stellar radius of 0.757 ± 0.003 Rq, and an inclination of 
85.61 ± 0.04 degrees. 

These data are clearly the best of their kind and we 
have attempted to fit them with our techniques and eq. 
[H The results are displayed in Fig. |8l The data are 
plotted as black hexagons, while the models are for var- 
ious values of P„. One model (dashed green) assumes 
10 X solar metallicity. All these models, save one, assume 
(Pq: Pi) = (0.1, 1.0) bars, not our default pair, but this 
makes little difference. As indicated in i j4.21 we can fit 
the contrast ratio in IRAC 4 at secondary eclipse rather 
easily, with a slight preference for a small non-zero Ko . As 
Fig. 121 suggests, a super-solar metallicity might also do 
the trick, but the metallicity dependence is rather weak. 
Data in other bandpasses should break the degeneracy. 

^ However, this is only the first of many anticipated. 



10 



However, we can not fit the small day/night difference 
with any of our models. The data seem to imply a se- 
vere degree of heat redistribution, one that is still not 
captured even with our P„ = 0.5 model. We note in 
passing that models with P„ — 0.5 do not imply that the 
dayside and nightside should look the same, only that 
the integral fluxes over the entire spectrum should be 
comparable. Since the dayside is irradiated, while the 
nightside emits into the blackness of space, as the upper 
right-hand panel of Fig. [1] indicates, the T/P profiles 
at a = 0° and a — 180° are different. This translates 
quite naturally into different day-night contrast ratio dif- 
ferences for different wavelengths, even for P„ = 0.5. 

What we seem to be seeing in the Knutson et al. 
(2007b) data are atmospheric inhomogeneities, thermal 
structures (vortices?), on the surface of HD 189733b. 
That the hot spot and the coolest spot are in the same 
hemisphere, separated by only ^--^45°, suggests our sym- 
metric models are inadequate to fit this phase curve. It is 
of paramount importance that a full light curve over all 
phase angles be taken in a variety of wavebands. Well- 
sampled data at longer wavelengths would be particu- 
larly welcome. Moreover, model phase curves need to 
be sophisticated enough to incorporate temporal and 3D 
spatial variations. The light curves and night-side heat- 
ing and thermal profiles depend centrally on jet streams, 
winds, and general thermal redistribution. This puts a 
premium on developing GCMs with reasonable global 
dynamics coupled to realistic radiative transfer models. 
Such models do not yet exist for the study of EGPs. 

6. DISCUSSION AND CONCLUSIONS 

In this paper, we have constructed atmosphere and 
spectral models for all the close-in extrasolar giant plan- 
ets for which direct-detection data from Spitzer have 
been published (except for the "Neptune" GJ 436b). 
These models incorporate the effects of external stel- 
lar irradiation, detailed atmospheres, heat redistribution, 
and, for some, a model for stratospheric heating. Com- 
paring the resulting suite of models with the data for 
these six EGPs, we have derived constraints on their at- 
mospheric properties. We find, as did Burrows et al. 
(2007b), that many severely irradiated EGPs can have 
thermal inversions at altitude which translate into qual- 
itative changes in 1) the planet/star contrast ratios at 
secondary eclipse, 2) their wavelength dependences, and 
3) day-night flux contrasts during a planetary orbit. Ab- 
sorption features can flip into emission features, plane- 
tary fluxes at long wavelengths can be enhanced, and the 
secondary-eclipse spectra in the near-IR can be altered 
significantly. What is more, we find a correlation be- 
tween the importance of such stratospheres and the flux 
at the substellar point on the planet. 

Hubeny, Burrows, & Sudarsky (2003) and Burrows, 
Sudarsky, & Hubeny (2006) showed that strongly irradi- 
ated atmospheres can experience a solution bifurcation to 
an atmosphere with an inversion for which water spectral 
features are reversed from troughs (absorption) to peaks 
(emission) . This possibility is supported by the good fits 
obtained by Burrows et al. (2007b) to the HD 209458b 
IRAC data (Knutson et al. 2007c) and by our models in 
this paper for the subset of irradiated EGPs for which the 
presence of stratospheres is suggested (in particular HD 
149026b and, perhaps, v And b). In Hubeny, Burrows, 



& Sudarsky (2003) and Burrows, Sudarsky, & Hubeny 
(2006), as well as in the prescient paper by Fortney et al. 
(2006), the absorber was gas-phase TiO/VO, which for 
hot atmospheres in chemical equilibrium can exist at low 
pressures at altitude and not just at high temperatures 
at depth. The upper-atmosphere absorber that is pro- 
ducing stratospheres for the higher values of Fp might 
indeed be TiO/VO, but a "cold-trap" effect can operate 
to deplete the upper atmosphere of TiO/VO. However, 
when mass loss is ongoing, as we know to be the case for 
HD 209458b (Vidal-Madjar et al. 2003,2004), the atmo- 
sphere is constantly being replenished and TiO/VO at 
some non-zero abundance remains a viable option. Such 
vigorous mass loss is expected for those planets with the 
highest values of Fp, and in this paper we have discov- 
ered a possible correlation between Fp and the existence 
of thermal inversions and stratospheres. Hence, the pos- 
sible mass- loss/stratosphere connection may make for a 
compelling scientific narrative. 

The tholins, polyacetylenes, or various non-equilibrium 
compounds discussed in the context of solar-system bod- 
ies could also be the necessary upper-atmosphere op- 
tical absorber. Given the stellar UV and integral-flux 
regimes experienced by strongly-irradiated EGPs, such 
species might be photolytically produced with sufficient 
abundance (Burrows et al. 2007b; Marley et al. 2007). 
However, to study these molecules requires a full non- 
equilibrium chemical network and we won't attempt this 
here. Clearly, what the high-altitude absorber actually 
is, TiO/VO or some other compounds, awaits investi- 
gation and is the primary reason we parametrized its 
opacity with Kc^. 

The trend with Fp we have uncovered suggests, how- 
ever crudely, that those EGPs with values of Fp higher 
than HD 209458b's (--10^ erg cm-^ s'^) may well have 
stratospheres. Table [T] provides the needed numbers. 
What this table suggests is that TrES-2, TrES-3, TrES-4, 
HAT-P-2b, HAT-P-4b, HAT-P-5b, HAT-P-6b, OGLE- 
TR-lOb, WASP-lb, WASP-3b, X0-3b, OGLE-TR-56b, 
0GLE-TR-211b, and OGLE-TR-132b, in addition to HD 
149026b, are strong candidates for having stratospheres, 
with all the consequences implied for their spectra and 
light curves (SjU [J5]). Close-in, though non-transiting, 
EGPs with high values of Fp (such as t Boo b, to name 
only one of many) are also hkely to have thermal inver- 
sions. They too should manifest the spectral discrimi- 
nants identified in Figs. |4] and [5] and the changes in the 
phase curves discussed in fj5] and suggested by Figs. [6] 
and [71 Note that by including in the above list HAT- 
P-2b, which has ~10x the mass of the average close-in 
EGP, we have not addressed the possible role of grav- 
ity in these systematics. Though we suspect gravity is 
sub-dominant when compared with Fp, the dependence 
of upper-atmosphere physics and chemistry upon gravity 
should prove worth exploring. 

It is unlikely that we will soon obtain spectral data for 
the OGLE planets. However, it is distinctly possible that 
EGPs listed in Table [T] with FpS slightly lower than HD 
208459b's will have "weak" stratospheres, as we specu- 

^ However, for some of the models in Burrows et al. (2007b) 
we used equilibrium TiO/VO abundances and the corresponding 
molecular opacities (Sharp & Burrows 2007). These models repro- 
duced the HD 209458b IRAC data reasonably well. 
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lated may be the case for HD 189733b. This could in- 
clude X0-2b, HAT-P-lb, WASP-2b, and, perhaps, XO- 
Ib. We note that Fp for TrES-1 is lower still (^0.43x10^ 
erg cm^^ s^^) and that this planet shows good evidence 
for water in absorption and no appreciable stratosphere 
( §4.31 and Burrows, Hubeny, & Sudarsky 2005). Hence, 
we have a hint at a breakpoint between EGPs with and 
without significant stratospheres and thermal inversions. 

Despite speculation to the contrary, our models with 
abundant atmospheric water are fully consistent with all 
the Spitzer data for all the measured EGPs, though at 
times the water is in emission, not absorption. This con- 
clusion is consistent with the possible identification of 
water in HD 189733b by Tinetti et al. (2007) (though 
see Ehrenreich et al. 2007) and in HD 209458b by Bar- 
man (2007). 

We find that the family of close-in EGPs probably 
boasts a range of values of P„ from ~0.1 to ~0.4. How- 
ever, our constraints on this parameter are rather weak, 
particularly given the possible complicating effects in 
some EGP atmospheres of a stratospheric absorber. Cur- 
rently, the presence of such an absorber makes it easier 
for values of P„ that are not small to explain the data, 
even for v And b. However, the magnitude of these ef- 
fects is hard to pin down with rigor, other than to say 
they are in evidence — there remains a slight degeneracy 
between P„ and "kc-" Without a more first-principles 
theory concerning the chemistry, spectroscopy, and abun- 
dance of this extra absorber at low pressures, and con- 
cerning the stellar environment in which it arises, use- 
fully constraining P„ may continue to be difficult. This 
is the case even if TiO and VO fit the bill, since their 
steady-state abundances would still be an issue. More- 
over, much better models for redistribution are urgently 
needed. Lunine & Lorenz (2002) have speculated that 
planetary atmospheres with jet streams to redistribute 



dayside heat to the nightside adjust their wind dynam- 
ics to maximize the rate of entropy generation. Their 
formalism suggests a value of P„ of ~0.2, which is not 
inconsistent with any of the currently known data on 
EGP atmospheres. 

The direct-detection data we have addressed in this pa- 
per are clearly only the first of many anticipated Spitzer 
contrast and light-curve measurements of strongly- 
irradiated EGPs. MOST will continue its campaign, 
future dedicated space missions will be proposed, and 
ground-based IR telescopes may have a role. Transiting 
EGPs are continuing to be discovered at an impressive 
rate that will not soon abate, providing an expanding 
catalog for follow-up and characterization. JWST is in 
the wings to revolutionize the field, and will come on line 
in the middle of the next decade. This paper is meant to 
provide a broad theoretical context for these initiatives 
and an interpretive vocabulary with which to address the 
ongoing study of extrasolar planets in tight orbits around 
their illuminating stars. 
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APPENDIX 

AN IMPROVED TREATMENT OF THE REDISTRIBUTION OF STELLAR IRRADIATION FROM THE DAYSIDE TO THE 

NIGHTSIDE 

In this Appendix, we describe a slightly more physical version of our previous treatment (Burrows, Sudarsky, & 
Hubeny 2006) of the day to night side redistribution. The total incident stellar flux (expressed as the iJ-moment) at 
the planetary surface is 

1 /RA^ cr 1 /R^ 



Hpxt — 



47r 



rp4 _ 

■ -'cff — 



(Al) 



where i?* is the stellar radius, T^s is the effective temperature of the stellar surface, and a is the planet-star distance. 
The basic feature of our model is an assumption that out of this total incident stellar flux, a fraction P„ contributes 
an additional source of energy on the nightside, and is removed from the dayside. P„ is bounded between 0.0 (no 
redistribution) and 0.5 ("complete" redistribution) and is conceptually the same as the redistribution parameter 
employed in Burrows, Sudarsky, & Hubeny (2006), but is here implemented slightly differently. We introduce 



2 



(A2) 



Formally, we take P„ > at the nightside (signifying a gain in energy), and P„ < at the dayside (signifying a sink 
of energy) . 

We define a local gain/sink of energy, D(m), such that 



D{m) dm — H^^ 



(A3) 



We assume that D{m) is non-zero only between column masses mp and mi (specified through limiting pressures Pq 
and Pi, which for our default set of calculations are 0.05 and 0.5, bars, respectively). We have studied the model 
dependence on the values of Pq and Pi around 0.1 to 1.0 bars. The results suggest only a modest effect on the overall 



12 



spectra at secondary eclipse. However, as the top-left panels of Figs. [T] and |4] suggest, one can optimally fit the IRAC 
2 to IRAC 1 flux ratio for HD 209458b by strategically placing the redistribution band between the corresponding 
photospheres, and thereby cool the IRAC 1 "r = 2/3" surface relative to the IRAC 2 "r = 2/3" surface. While this did 
not motivate our default values of Pq and Pi (used, one notes, for all objects in this study, not just HD 209458b, and 
motivated by the desire to redistribute heat near the TRossciand = 2/3 level), the reader should be aware that detailed 
fits to the secondary eclipse data depend on the choice of Pq and Pi. Clearly, better models of heat redistribution 
than we have employed here are called for. 

With this caveat in mind, we consider two models for D{ni). The first one assumes a constant D{m) between mg 
and mi. 

D{m) = , (model 1) . (A4) 

mi — mo 

The second model, model 2, assumes a linearly decreasing D{m) between mg and mi, in such a way that D(m) reaches 
at m = mi. This is our default model for the calculations of this paper and its functional form is: 

„ , N 2iJirr mi — m / , , s , . , 

D{m) = , (model 2) . (A5) 

mi — Too TOi — Too 

Therefore, D[m) is non-negative on the nightside and is non-positive on the dayside. 
The first two moments of the transport equation read 

^ =«:,(J,-B,), (A6) 
am 

and 

^ = X.H,, (A7) 
dm 

where k^, is the absorption coefficient per gram, is the total extinction coefficient (absorption + scattering) , and m 
is the column mass, related to the pressure by the relation P = mg, where g is the gravity. 

Since we stipulate sinks or sources of energy at certain layers, the usual radiative equilibrium (or radia- 
tive-|-convective) equilibrium does not apply. Instead, it is replaced by the following energy equation, which can 
be written in two different ways: 

a) using the frequency-integrated first moment of the transfer equation. 



/•OO 

/ Hv{Jv — Bi,) dv = —D{m) 
Jo 



(A8) 



because the energy gained per unit mass, D(m), is balanced by the net radiation loss per unit mass, given by the 
integral on the left-hand-side of eq. (|A8|) . 
b) using the equation for the frequency-integrated flux, H = H^dv. From eqs. ()A6P and (jA8|, we have 

^ = -i.(TO), (A9) 



which we rewrite as an equation for the integrated H as 



H{m)^H^s+ D{m')dm' , (AlO) 

where to^ is a sufficiently large column mass at which one has D{md) = (that is, deeper than the region of the 
sources/sinks), and 

^eff ^ ^Tltt (All) 

is the nominal total flux deep in the atmosphere, expressed through an effective temperature, Ti„f 

In view of a simple linear form of D{m), the integral in cq. (jAlOp can be evaluated analytically. We obtain for 
model 1: 

H{m) = HcS , m > mi 

H{m) ~ Hcs + -ffirr , rno < m < mi 

TOi — Too 

H{m) = Hcff + Hirr m < mo , (A12) 



and for model 2: 



H(m) = Hcs , TO > TOi 

H(m) = Hctf + Hirr i — I , 1TlO ^ '77 < mi 

V mi - mo/ 

H{m) = Hcff + Hi„ m < Too . (A13) 
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As suggested by Hubeny & Lanz (19 95), i t is numerically advantageous to represent the energy balance equation as a 
linear combination of eqs. (jASP and (jAlOp . where the ff- moment is expressed as Hi, — dK^/dm — d{f^J^)/dm, where 
fy is the Eddington factor. Using the Eddington factor enables us to consider only one radiation moment, Ju, as an 
unknown quantity. The Eddington factor is not taken as an unknown; instead it is held fixed in linearization and is 
recalculated in the formal solution step. 

The above equations do not consider convection, so they apply in the radiative zone. In the convection zone, eqs. 
(jA10[) . (|A12p . and (jAlSp remain unchanged, provided we take 

^conv 

H.^ = - ^ , (AM) 

where is the convective flux. Equation (|A8p is modified in the presence of convection to read 



1 dF'^°'^^ 

K,[J,~B,)dv^-D{m)-- . (A15) 

47r dm 

CoolTLUSTY solves the energy equation specified through eqs. (jAlip - (jAlSp numerically and self-consistently 
with the set of equations of radiative transfer. 

However, one can gain deeper physical insight by developing a simplified gray model, in which we can actually derive 
analytic expressions for the local temperature. More importantly, such a model allows us to understand what values of 
the limiting pressures for the sink region on the dayside and the corresponding optical depths are physically acceptable. 
We develop these expressions in Appendices iiB]& below. 

A SEMI- ANALYTIC, GRAY MODEL WITH REDISTRIBUTION 

In CoolTLUSTY, the exact energy balance equation is solved self-consistently with the radiation transport equa- 
tion. However, it is very useful to develop a simple gray model that allows us to study the conditions under which 
the structural equations have a solution at all. On the nightside, there is always a solution because we are adding 
energy at certain layers. However, on the dayside, we remove energy at certain layers. If we require that these layers 
are too deep in the optically-thick part of the atmosphere, the only way energy can be removed is to create a negative 
temperature gradient. If, moreover, we require the region with a negative gradient (the region of the energy sink) to 
continue to depth, we would eventually reach a negative temperature, which is clearly unphysical. 

To demonstrate this and to construct an analytic atmosphere model with redistribution, we closely follow the 
derivation of the analytic model given in Hubeny, Burrows, & Sudarsky (2003), generalizing it to account for departures 
from radiative equilibrium due to the stipulated sources/sinks of energy. We write the frequency-integrated moment 
equations (jA6p and (|A7p using the mean opacities: 

dM 

— ^njJ-KBB, (Bl) 
dm 

and 

dK , , 

where kj, kb, and xh are the absorption-mean, Planck-mean, and flux-mean opacities, respectively, and J,H,K are 
the frequency-integrated moments. 
The energy balance equation reads (neglecting convection) 

KjJ - kbB = -D{m) , (B3) 

which is just another form of eq. (|A8[) . 
First, we obtain the solution for the second moment K. We write the second-moment equation as 

dK(T) , , , , 

= H{t) , (B4) 

where dr = XHdm « Xmssdm is the flux-mean optical depth, which can be approximated as the Rosseland optical 
depth. In order to integrate the second-moment equation analytically, we introduce the limiting optical depths tq = 
r(rno) Stiid ti = r(TOi), and adopt the following approximation for H{t): 

H{T)^H,s + Hi,J^^^) , forTo<r<ri, (B5) 

together with the exact expressions H{t) = H^ff for r > ti and H{t) — H^s + Hi„ for r < Tq. Here, n = 1 for model 
1, and n = 2 for model 2. This equation is easily solved and yields: 

K{t) = K{0) + (H.s + H,rr) r, r < tq 



K{t)^K{0) + H,bt + H,,Ato 



Tl - Tq 



1 



Tq <T <Tl 



K{t) = K{0) + H,sT + Hi,, [ ^— Tl + —— To] , r > n . (B6) 
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With the LTE-gray model, this equation is in fact an equation for the local temperature. We invoke the Eddington 
approximation, K = J/3, and use the energy balance equation 

kbB = KjJ + D{m) , (B7) 

and the fact that B — (cr/7r)r^ to derive the local temperature. We introduce the following quantities: 

J(0) = 3i^(0) = aJoxt , i^Gxt = fuJe^t , -ffirr = Pn^ext, and W = iJoff/^fext , (B8) 

and we write D{t) instead of D{m) in eq. (jB7[) as an approximate expression 



D{r) 



(B9) 



Tl - To Kb Tl - To 

where x and kb are the average values of the flux- mean and the Planck- mean opacities in the interval (To,ri). Note 
that in the strict gray model, e = 1. Note also that the ratio w is given by 

w = (ri„t/r,)4 2 (i?,/a)-2 , (BIO) 

so that for the case of strong irradiation, it; <C 1. Since here we are interested in deep layers, we use the approximation 
ftj — Kb, (which is, however, not valid at the surface layers). Using the above defined quantities, we can express the 
integrated Planck function (i.e., temperature) through H^xt only and derive: 

B= f3 + wT + 3P„g(T) + i/ext, (Bll) 

Tl - To _ 

where we denoted (3 = a/{3fH) (which is a constant of order unity), and 

1 



9(t) = To + 



1 



(ti - To) 



1 - 



Tl 



Tl - To 



n+1 



(B12) 



for To < T < n; q{T) — for r < tq, and (?(t) = q{Ti) for r > ti. Notice also that since B = {a/n)T'^, eq. (|Bll[) can 
be understood as an equation for the local temperature. 

CONDITION FOR THE EXISTENCE OF THE SOLUTION ON THE DAYSIDE 

We now investigate the existence conditions for the solution on the dayside. We first make the following approxi- 
mations: a = 2, fn = 1/2; thus, /3 — 4/3, e = 1. We assume that tq ^ ti, so we neglect it in the expression for ^(t). 
We also introduce the notation p„ = — Pn, noting that p„ is a positive quantity. Let us first take model 1, in which 
n = 1. Using all the above approximations, we obtain for the region between tq and ti: 



B{t)/H,-_ 



Pn 
Tl 



+ WT - p„r 1 



2ti 



(CI) 



Since t\ is typically larger than 1, we neglect the term p„/ti compared to /3. This is not necessary for the formal 
development, but it simplifies the resulting expressions. The derivative dB /dr — w — p„(l — t/ti); therefore, the local 
minimum of B{t) is at t = r^in = ti(1 — w/pn). Since the most interesting case is for strong irradiation where w 1, 
we see that the minimum of B is close to ri. The value of B{t) at the local minimum is 

Bmin = i?(T„,i„) = /3 - ^ f 1 - . (C2) 

2 V PnJ 

The condition for the existence of the solution is that Smin > 0, and, thus. 



2/3 / w 

Tl < 1 

Pn \ Pn 

Taking the most interesting case of strong irradiation, we can neglect the second factor, and write simply 

2/3 8 

Tl < — « . 

Pn iPn 

An analogous analysis for model 2 gives a similar condition, 

4 

Tl < — . 

Pn 

For instance, this demonstrates that for pn = 0.5, the case with the maximum energy sink on the dayside and the 
maximum degree of redistribution, the deeper (high-pressure) limit of the sink region must be at optical depths less 
than 8. Eq. IC5I indicates that this limit is larger for smaller j>„ (e.g., it is 40 for p„ — 0.1). Hence, we have derived 
consistency conditions for our redistribution algorithm that have physical content and in our choices for Pq and Pi , we 
are careful not to exceed this condition. Our default values of Po and Pi, 0.05 and 0.5 bars, respectively, translate into 
an optical depth range of a fewxO.l to '^a few, within the consistency constraints for all the P„s. Indeed, if we were to 
exceed the consistency constraints in our atmosphere calculations, they would not converge and the simulations would 
crash numerically. 



(C3) 



(C4) 



(C5) 
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THE ORIGIN OF THE / = 2/3 TERM 



To derive the proper / factor, we now expand upon the formahsni of Appendix fJX] The total energy flux received 
by a unit area on the planetary surface at angle Oq from the substellar point is given by 

F{^io) = ^TT (^^^ H.fio, (Dl) 

where = cos 60 and 9q is the angle between the normal to the planetary surface and the direction toward the star. 
For simplicity, we assume that the angular diameter of the star is small. Consequently, all rays coming from the star 
are parallel. iJ* = -F'*/47r and is the radiation flux at the surface of the star. The average flux received by the 
planet is then 

i^av = ^ F{no) dfio = I ^ttH, (^^^ , (D2) 

which explains the origin of the f ~ 1/2 ansatz. 

To improve upon this, we assume the Eddington approximation, J = 3K, and ignore convection, which allows us to 
write down an analytic solution of eq. (B4) for J: 

J{t) = Jo + 3iJr , (D3) 
where AnH is the interior planetary flux. To obtain the constant Jq = J(0), we employ the formal solution of the 
transfer equation for the specific intensity: 

/•oo 

/(0,m)=/ {Kj/KB)J{t)e-'/^dt/fi, (D4) 



since the source function, S, is given by 5 = i? = / hb)J ■ As shown by Hubeny et el. (2003), Kj/nB can differ 
significantly from unity, but only for low values of the flux- mean optical depth, r ^ 1. Therefore, we set kj/kb — 1, 
and, using eq. (|D3p . we integrate eq. (|D4p to obtain 

7(0, m) = Jo + 3i7/^ , (D5) 
which is the well-known Eddington-Barbier relation. The mean intensity on the planetary surface is given by 

J(0) = (1/2) 7(0, m) d^x + (1/2) y" P^\^l)d^l , (D6) 

where 

P^\^Ji) = 5{-^i-^^o)F{^^^)/^, (D7) 

because we have assumed all incident rays from the star are parallel. 5{) is the Dirac J-function. Therefore, 

Jo = J(0) = ijo + ^i7+^, (D8) 

and, consequently 

j„^3i^+^M. (D9) 

TT 

The specific intensity is then given by 

/(/i,Mo) = ^^+3i7(M+l), (DIO) 

TT 

where we have given the explicit dependence of the emergent specific intensity on /iq- 

For close-in planets, the irradiation flux is much larger than the intrinsic flux, AttH , so we neglect the second term. 
The local atmosphere characterized by angle ^0 exhibits, within the present approximations, an essentially isotropic 
emergent radiation pattern, independent of the local polar angle n and dependent only on the angular distance from 
the substellar point, ^o- For the total planetary flux close to secondary eclipse received by a observer at a distance D, 
we have the expression 

' If',., 4/i?,^' 



{D/RpyFohs = / /(Aio,Mo)Moc^Mo ~ J ^(Mo)MocfAio = ^ y— ) ^* ■ (DH) 



2 



For an average atmosphere characterized by the parameter, /, the external flux is given by eq. (jD2[) . where we replace 
1/2 by /. We obtain 

and, therefore. 



/(Ai,Mo) = ^^=4(^) 77,/, (D12) 



(7?/Rp)2Fobs = y I{iJLo,ti(,)^io duo = 4 77, / y /xo d^o - y ■ (D13) 

In order to get agreement between eqs. (|D11|) and (jD13|) . we have to set / = 2/3. This is the origin of our use of this 
value. 
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TABLE 3 

Spec^tral. Photometric, and Composition Measurements of Extrasolar Giant Planets 
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HD 180733 b 



7.5 14.7 
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Spitzor/IRS 

Spitzer/IRAC4 
Spitzer/IRS 

Spitzcr/IRACl 

Spitzcr/IRACa 

Spitzci7IRAC4 

Spitzcr/IRAC2 

Spitzcr/IRAC4 

HST/STIS 

HST/STIS 

HST/STIS 

HST/STIS 
HST/STIS 
MOST 

MOST 
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Spitzer/IRS 



ad. 



Second. Eclipse 

Spoctrum 

0.003392(55) 
0.00551(30) 

0.02356(020) 

0.02436(020) 

0.0239(02) 

0.00066(13) 

0.00225(36) 



< 1.34 X 10-"* (,3o-) 

< 4.88 X 10-" (Icr) 

< 3.9 X IQ-s (la) 

< 0.0003 (la) 



0.00094(9) 
0.00213(15) 
0.00301(43) 
0.00240(26) 

Spectrum 



AFp/F, i , 



CmrirTipntn 



No H2O7, CH^l 



Grillmair ot al. (3007)^ 



TrES-1 

HD 209458 b 



HD 149026 b 
V And b 
HD 179949 b 
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51 Peg b 



8.0 
16.0 

3.6 

5.8 

8.0 

4.5 

8.0 

Ly a 

0.12-0.17 

0.3-0.5 

0.58-0.64 

~0.94 

0.4-0.7 

0.4-0.7 
2.2 
3.6 
8.0 

3.6 
4.5 
5.8 
8.0 

7.5-13.2 



8.2-13.2 


Spitzer/IRS 


Spectrum 


24.0 


Spitzer/MIPS 


0.00260(46) 
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Spitzer/MIPS 


0.0033(3) (?) 


8.0 
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0.00084(11) 
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0.0012(2) 



< 0.015 (2a) 

< 0.0015 {2a) 



0.0029(7) 
< 0.019 {2a) 
0.00141(33) 



< 0.017 {2a) 

< 0.0007 (2a) 



Light curve 
Peak-up mode 
Primary Transit Depth 
(Tinetti et al. 2007) 

H2O identified — 
(Burrows et al. 2005) 
HI, Rp = 4.3 Rj 
C, O 

Rp = 1.3300(6) Rj 
at Balmer cont.? 

Na D (transit) 
H2O identified (transit) 

As 
A, 



< 0.68 (3a) 
■ < 0.25 (Icr) 
~ 4.0 ± 4.0% 



Light curve 
Light curve 
P„ > 0.32 (la) 
Temperature Inversion - 
H2O in emission — 
(Burrows et al. (2007b) 

Features at 
7.78, 9.65 um ? 
No H2O, CHi? 
Fp = 0.40 ± 0.19 
mjy at 12 /im 



Ligfit curve 
Light curve 
Light curve 
Pn < 0.30 (la) 



Light curve 
Light curve 



Knutson et al. (2007b) ^ 
Deming et al. (2006) 
Beaulieu et al. (2007) 
Beaulieu et al. (2007) 
Knutson et al. (2007b) 
Charbonncau et al. (2005) 
Charbonncau et al. (2005) 
Vidal-Madjar et al. (2003)"* 
Vidal-Madjar et al. (2004)5 
Ballester et al. (2007)"^ 

Charbonncau et al. (2002)'' 
Barman 2007 

Rowc et al. (2006)^ 
Rowc ct al. (2006) 
Rowc ct al. (2007) 
Richardson ct al. (2003) 
Cowan et al. (2007) 
Cowan et al. (2007) 

Knutson et al. (2007c) 
Knutson et al. (2007c) 
Knutson et al. (2007c) 
Knutson et al. (2007c) 
Richardson et al. (2007)^ 



Swain et al. (2007)1° 

Deming et al. (2005) 

Deming (private communication) 

Harrington et al. ^2007) 

Harrington ct al. (2006) 

Cowan et al. (2007) 

Cowan et al. (2007)^ 

Deming ct al. (2007) 

Demory et al. (2007) 

Cowan et al. (2007) 

Cowan et al. (2007) 



OSI)LE-4]i3ofe= -2r2is tabic c 

observations, the planet to star 
and the reference. Numbers in h 
data (except the 8-/im point for 
argue that a flat spectrum is ea 
eelip.so data 



^MOiT/^jriLwing eoQuQai:7(S^ 4&0-^Zi the pi 
lux ratio during the secondary eclipse (Fp/F*), t] 
old are measurements (not upper limits) at secon 
GJ 436b at secondary eclipse) arc presented in this pap 
ised by the lack of significant water or methane absorpt 
>f Beaulieu et al. (2007) and Knutson et al. (2007b); 



t, spectral region of the observations, the telcsc' ope \Snellenir&rCfi)SfitlG)s((<£@tG7t)he 
:Lmplitude or peak to trough variations of the light curve (AFp/F*), brief comments, 
y eclipse or actual light curve (AFp/F+) measurements. Models for all these bolded 
paper. Numbers in parentheses are errors in the last digits. rillmair et al. (2007) 



an. However, this is not consistent with the cor 
nutson et al, (2007b) also derived the longiti 
Lie the coolest region was shifted about 30 degr. 
asit radius at 8 fim of 1.137 ± 0.006 R j . a stella 



elusions of Tinetti et al. (2007), 
idinal dependence of the surface 
icis west of the anti-stellar point, 
r radius of 0.757 ± 0.003 R^-. 



brightness and found a hot spot shifted by 16 zb 6 degrees east of the substellar poir 
They also found an indication of nonzero eccentricity with ecos = 0.0010 ± 0.0002, 
an inclination of 85.61 ± 0.04 degrees, where uj is the longitude of periastron; owan et al. (2007) constrain P„ for HD 209458b and HD 179949b to be P^^ > 0.32 (Icr) 
and < 0.30 (Ict), respectively; Pn is the fraction of the total energy incident on the dayside of the planet which is transfered to and radiated out on the night side of 
the planet; idal-Madjar et al. (2003) detected atomic hydrogen in the planet's atmosphere with a transit absorption depth of 15 ± 4(lo-)%, and evaporation at the rate of 
> 10^*^9 idal-Madjar et al. (2004) detected oxygen and carbon in the planet's atmosphere; allester et al. (2007) claim to have identified HI absorption in the Balmer 

continuum, but Barman (2007) challenges this interpretation; harbonneau et al. (2002) found that the transit depth at the Na D feature is deeper by 2.32 ± 0.57 X 10~'* 
than in the continuum, which is interpreted as a detection of Na in the planet's atmosphere; owe et al. (2006) constrained Ag < 0.25 (la-), or Ag < 0.68 (3cr), but Rowe 
et al. (2007) constrained Ag to be below 8% to 1 — ct. Ag is the geometric albedo in the optical; ichardson et al. (2007) claim to have detected a broad emission peak 
centered near 9.65 /L*m which they attribute to the emission by silicate clouds, and a narrow unidentified emission feature at 7.78 ^^m. They say that models with water 
absorption fit the data poorly. However, Burrows et al. (2007b) conclude that water is in fact seen in emission, not absorption; wain et al. (2007) determined the planet 
flux at 12 jum to be 0.40 ± 0.19mjy and the normalized secondary eclipse depth to be 0.0046 ± 0.0006. They are updating their absolute calibrations. 
1 „2 j,3 ^4 „5 „6 „7 ^8 „9 „10 „ 
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Fig. 1. — Temperature-pressure profiles for the six close-in planets studied in this paper. Dayside profiles incorporate the external 
substellar irradiation/fiux given in Table [T] and an internal flux for the planet corresponding to the temperature of 75 K at the lower 
boundary. A sink of energy corresponding to the particular value of P„ was introduced between pressures of 0.05 and 0.5 bars. The 
nightside is calculated without irradiation, assuming an energy source corresponding to the same value of P„ at the same pressures 
employed for the dayside sink. Entropies at the bases of the convection zones on the day- and night-sides were approximately matched by 
adjusting the internal planetary fiux at the bottom of the nightside model at the same gravity. Models with an extra upper-atmosphere 
absorber in the optical are included for HD 209458b, HD 189733b, HD 149026b, v And b, and HD 179949b. See text in f||]for a discussion 
of these panels. 
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Fig. 2. — "Brightness" or "formation" temperature spectra on the dayside for three models (Pn = 0.1, 0.3, and 0.5) of TrES-1. The 
formation temperature at the particular wavelength is the temperature where the optical depth at this wavelength reaches 2/3. See text 
for a discussion. 
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Fig. 3. — Spectral energy distribution (SED) {XFx versus logio(A)) for the dayside of TrES-1 for three values of the day-night heat 
redistribution parameter P„. Notice that observations in IRAC and MIPS bands cover only a small fraction of the SED of the planet. 
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TrES-l 

(At Secondary Eclipse) 
Charbonneau et al. (2005) 




HD149026b 

(At Secondary Eclipse) 
Harrington et al. (2007) 



P„=0.0, K,=0,00 cm' 

P^ = 0.0, «:^=0.6 i 
P„=0.3, (C =0.00 ' 




Fig. 4. — The planet /star flux ratios versus wavelength from ~1.5 fim to 30 fim for various models of four transiting EGPs measured by 
Spitzer at secondary eclipse. Notice the different scales employed in each panel. Models for different values of P„ (EJ and fte are provided 
where appropriate and the data from Table \3\ for each planet are superposed. The plot legends indicate the color schemes used for the 
different EGPs. On the upper- left panel (HD 209458b), models with the lighter gray shade(s) are for the higher value(s) of P,i. Notice 
also that two different values for the flux at 24 fira (green) are shown on this same panel. The one with the question mark is a tentative 
update to the Deming et al. (2005) 24-fim measurement, kindly provided by Drake Deming (private communication). If the flux at 24 fim 
is indeed ~0.0033±0.0003, then our model(s) with inversions provide the best fit at that wavelength as well. Note that the comparison 
between model and data should be made after the band-averaged flux-density ratios of the detected electrons are calculated. This not only 
incorporates the significant widths of the 5piteer bandpasses, but the fact that one should compare photon counts (or detected electrons), 
and not monochromatic fluxes. The result is that the theoretical IRAC predictions do not actually vary on the figures as much as do the 
plotted spectra and that the predicted contrasts, for instance between IRAC 1 and IRAC 2, are more muted, even when the pronounced 
bump at and near ~3.6 /im obtains. However, to avoid the resultant clutter, we do not put these bandpass predictions on these spectral 
plots. See text for a discussion of each irradiated planet and the inferences drawn. 
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Fig. 5. — The same as Fig. |4l but for various models of HD 209458b between 1 fim and 4 fim. The dependence on both P„ and the 
presence and strength of a thermal inversion is greatest in this wavelength region. Note that a thermal inversion flips what would be water 
absorption features into emission features, altering the interpretation of any data in this spectral region significantly. Superposed is the 
Knutson et al. (2007c) data point for IRAC 1 and the Richardson et al. (2003) upper limit near 2.2 fim. See text for a discussion. 
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Fig. 6. — Theoretical light curves for the non-transiting EGP v And b in the MIPS 24-/xm, band for different inclinations [i = 45° and 80°), 
values of P„ (0.0 and 0.3), values of (0.0 and 0.2 cm^/g), and a range of planetary radii. Superposed are the light-curve measurements of 
Harrington et al. (2006). Note that Harrington et al. (2006) obtain only relative contrast values, not absolute values. Therefore, the data 
on each panel are shifted in absolute contrast space, while maintaining the measured relative values. The day/night contrast difference is 
preserved. In this way, we find the corresponding best fits consistent with the observations, but for different inclinations, etc. Note that 
the scales for the different panels can be different, though the right panel that faces each left panel has the same scale as that left panel. 
The left panels depict models without a thermal inversion, while the right panels depict models that have thermal inversions created using 
Ke = 0.2 cm^/g. Together, this set of panels illustrates the dependence on the three most important free parameters: P„, i, Hp, and fto. 
See text for a discussion. 
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Time (days) Time (days) 

Fig. 7. — The same as Fig. |6] but for HD 179949b in the IRAC 4 band for two different inclination angles (45° and 80°), two values of 
Pn (0.0 and 0.3), a range of planetary radii, and two values of Kc (0.0 and 0.08 cm^/g). The light curve data from Cowan et al. (2006) are 
superposed on each panel. See text for a discussion. 
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Fig. 8. — A comparison between the light-curve measurements of HD 189733b in the IRAC 4 band (8 /im) performed by Knutson et 
al. (2007b, hexagons) and our theoretical light curves for various values of Pn (0.1, 0.3, 0.5). Most of these models employ values for the 
redistribution pressure range (Pq and Pi) of 0.1 and 1.0 bars, but one model (dotted, and P„ = 0.3) uses (Pq, Pi) = (0.05, 0.5) bars. See 
figure legend for model parameters. Also included is a model with lOxsolar metallicity (dashed). All the P„ = 0.3 models are in green. 
See text for a discussion. 



